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Abstract 

We exploit a slightly noncollinear second-harmonic cross-correlation scheme to 
map the 3D space-time intensity distribution of an unknown complex-shaped ultra- 
short optical pulse. We show the capability of the technique to reconstruct both the 
amplitude and the phase of the field through the coherence of the nonlinear interac- 
tion down to a resolution of 10 fim in space and 200 fs in time. This implies that the 
concept of second-harmonic holography can be employed down to the sub-ps time 
scale, and used to discuss the features of the technique in terms of the reconstructed 
fields. 
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1 Introduction 



The study of ultrafast phenomena has been a major scientific priority during 
last decades covering different topics such as the study of radiation-matter 
interactions (1), transient response of molecules and atoms (JiJ), coherent con- 
trol of chemical reactions & or communication and inform^ ion technology 
(0). The growth of this field relies upon the development of sources of fem- 
tosecond radiation and of appropriate techniques able to provide time domain 
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information in the femtosecond scale. However, during the interaction of short 
optical pulses with a nonlinear medium, different mechanisms can lead to its 
reshaping into complex spatio-temporal structures with non-trivial light dis- 
tribution (Ja). As a consequence, their complete characterization requires a 
method capable of acquiring a snap-shot of their intensity distribution in the 
whole 3-dimensional (3D; x,y,t — z/c) comoving frame. 



Most of the available methods for pulse diagnostic provide information of the 
WP characteristics in a space of reduced dimensionality. The use of frequency 
resolved autocorrelation techniques (i.e. FROG, SPIDER) allows for the re- 
covery of the temporal intensity and phase profile of a given pulse but assumes 
uniform transverse spatial distribution (0;0). On the contrary, the character- 
ization of transversally localized beams often relies upon the optical imaging 
onto CCD cameras, therefore the temporal information is lost because of their 
integration times unavoidably larger than the optical pulse duration. Recently, 
a space-time characterization method based on extended SPIDER technique 
has been developed capable of resolving electric field characteristics in time 
and along one spatial coordinate (0). A quite direct way of obtaining spatio- 
temporal intensity profiles of a WP is to perform measurements with a streak 
camera, which allows a temporal resolution up to fractions of ps (|9|). This 
technique allowed the investigation of the dynamics of the breakup along the 
pulse envelope of a large elliptical beam propagating into a saturable Kerr 
nonlinear medium (flfl 11). However, also in this case, the space-time maps 
are intrinsically two dimensional (1 spatial + temporal dimensions). 



A different approach to the problem considers the retrieval of the pulse shape 
through an all-optical processing by means of spatially resolved detection sys- 
tems combined with gating techniques. The principle of the method is that of 
characterizing with spatial resolution an optical field that is proportional to 
the product E (x, t)E R (yi, t), where Eq{x, t) is the object to be measured and 
En(x.,t) is a suitable reference pulse. Since the product is different from zero 
only on the intersection of the support of both fields 1 , by translating the refer- 
ence with respect to the object, we get the possibility of recording information 
from different parts of the object. Among the linear time gating techniques, 
light-in-flight holographic recording has been the first technique which per- 
mitted the recording of dynamically evolving light fields during propagation 
(HI; IH 14). Recently, this technique has been adapted to study the propaga- 
tion of a 3 ps long pulse in linear media |l5h. Linear probing techniques were 
also exploited in order to obtain time resolved imaging, like the probing of the 
birefringence properties of plasma by means of a delayed, spatially-extended 



Rigorously, the support of a function is the set of points where its value is different 
from zero. For realistic optical fields that have exponentially decaying tails we can 
define a "practical" support, defined as the set of points in the space-time in which 
the field amplitude is larger than 1/e times the peak value. 
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100 fs pulses to investigate the dynamics of laser pulse focusing in air (Jl 

Nonlinear processes have been employed since long ago to resolve in time the 
evolution of ultrafast phenomena. Among them, the quadratic nonlinearity 
has been proved to be particularly versatile due to the fact that it provides 
easily terms containing the product of two optical fields. Recently, a type II 
degenerate parametric amplification scheme has been employed to obtain time 
resolved 2D images of a ps-pulse hitting a diffusing screen with 35 ps resolution 



(1171 ). thus yielding a 3D imaging. The same technique was later used to image 
an object embedded in a thick diffusing sample (1181). Althought our setup is 
actually an improved version of that described in §18), we point out that our 
conceptual approach is different from the study of the propagation of a wave 
front. In fact, in our case the propagation variable is fixed. 

Our goal in this article is to demonstrate the potentiality of the optical gating 
technique to acquire a high resolution space-time map of short, focused WPs 
in their comoving reference frame. Furthermore, we devise the capability of 
the technique to reconstruct both the amplitude and the phase of the WP 
thanks to the coherence of the nonlinear interaction. We propose a method 
that is based on quadratic type I interaction in a sum-frequency generation 
scheme either by a non-collinear second harmonic generation or by a collinear 
sum-frequency scheme. The latter has been used in (jlih . Here we discuss the 
first option, showing that if the interaction angle between the two interacting 
fields is small enough, then a reliable space-time map of the object pulse can be 
obtained. This can be achieved if the duration of the gate is much smaller than 
that of the object to be imaged. A holographic interpretation of the method 
permits to gain insight into the process of up-conversion of the space-time 
slices of the object into the SF field, and to prove that the coherence of the 
SF process is able to reconstruct the wavefront in both amplitude and phase. 
Our results confirm this possibility. The theoretical discussion of the method 
is followed by section 3, where we present the set-up and the experimentally 
reconstructed space-time intensity profiles of a parametric spatial soliton ex- 
cited by a 1 ps light pulse. For our setting, we estimate a mapping resolution 
of 200 fs in time and about 10/xm in space. The features of the technique are 
presented in section 4, pointing out the limitations that may arise and dis- 
cussing the possible implementations in each case. In the last section the main 
conclusions are presented. 



2 Description of the technique: intensity and field reconstruction 



In this section we explicitly show how a non-collinear sum-frequency (SF) 
scheme can be exploited to get high resolution space-time intensity maps of 
an unknown light wave packet with a space-time structure (object wave). The 
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recovery of a 3D intensity map is obtained by means of a short reference pulse 
which provides a time gating inside a nonlinear (NL) crystal, and generates 
a SF signal containing the information about a set of 2D slices of the object 
obtained by changing the reference delay. We first discuss the case for the re- 
construction of the object intensity profile, and then we show how the intrinsic 
coherence of the SF process allows in fact for a truly holographic recording of 
the unknown object. 

2. 1 3D Intensity profile mapping 

Let us denote the object (E ) and reference (E R ) wave packets as follows: 



where the complex functions Eo{x,y,t,z) and Ept(x,y,t,z) are the slowly 
varying envelopes of two waves with frequencies u± and cj 2 - Note that in this 
form the equations describe two wavepackets propagating in the positive z 
direction and colliding at an angle 9 = 2 arctan(/c x //c z ) in the x — z plane 



(here k — yk% + k% — 2ir/\ ). For a SF generation process occurring inside a 
quadratic nonlinear crystal, the polarization source giving rise to the SF can 
be written as: 



where the phase and energy matching conditions k z {ujz) = k z {uj\) + k z (u2), 
k x {oJ\) = —k x {uj2) and co>3 = uj\ + u 2 have been used. The SF field propagates 
along z direction and has a slowly varying envelope function that we will indi- 
cate by E S F(x,y,t, z). Now we introduce the following assumptions: 1) small 
depletion of both the Eq and E R fields; 2) negligible diffraction and dispersion 
within the crystal; 3) equal group velocities of the object, reference and sum- 
frequency pulses, namely uo, ur and usf, that is uo = ur = usf = u. Note 
that all these assumptions approximatively hold as long as the thickness of the 
crystal Az is small compared to the characteristic lenghts of the system (non- 
linear length, dispersion and diffraction lengths, pulse walkoff length). These 
assumptions allow to find a travelling reference frame for all the propagating 
pulses by introducing the retarded time r = t — z/u. The general partial differ- 
ential equations describing the interaction process then reduces to an ordinay 
differential equation for the envelope E SF . If we also introduce a time delay 
Tj on the reference wavefront, the equation takes the form: 



E = E (x, y, z, t) e *[«i*-M«i)*-M«i)*] + c . c . 
E R = E R (x, y, z, f) e *M-M<*)*+M«*0*] + c . c . 




(2) 




P SF oc 2E E R e i[ ^ t - kAu)3)z] + c.c. 



(3) 



dE SF (x,y,T,z) 
dz 



i2aE Q (x, y, r, z)E R (x, y,r - n, z) 



(4) 
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where a is the nonlinear coupling term. The last equation is readily integrated 
and, if the mixing crystal is placed at position z , it reads: 



E SF (x, y, r, z Q ) = ia2E Q (x, y, r, z )E R (x, y,r - n, z )Az 



(5) 



A deeper discussion about the meaning of this expression will be given in the 
following subsection. Here we just point out how the intensity profile of the 
SF field can be used to retrieve the 3D intensity map of the object. More 
precisely, for a given lag time Tj, the spatially dependent SF fluence profile 
(the CCD signal) S(x,y,Ti, zq) recorded just at the exit face of the mixing 
crystal is given by: 



S(x,y,Ti,z ) ~ (aAz) 2 / \E Q (x, y, r, z )\ 2 \E R (x, y, r - n, z )\ 2 dr (6) 



This expression provides the convolution between the intensity profiles of the 
object and reference wave packets, lagged in time by r«, and shows that the 
signal recorded by a CCD sensor is a linear function of the intensity of both the 
object and reference wavepackets. In the particular case in which the reference 
wave is spatially homogeneous in the transverse x — y plane, and temporally 
much shorter than the object, we can write expression 6 in the following form: 



where I Q = \E \ 2 - 

By imposing a set of delays (r,, i = l...n) to the reference WP with respect to 
the object, a reliable 3D reconstruction of the WP structure can be achieved 
by the collection of the n images S(x, y, 7$, z ). By changing the plane z , the 
temporal evolution of the WP can also be obtained. 

Notice that the model described above does not take into account the disper- 
sion of the mixing crystal. Therefore, the model as it is predicts no limits for 
the resolution of the map as long as arbitrary short reference pulses are avail- 
able. Actually the real mixing crystal has a finite bandwidth that limits the 
spatiotemporal resolution of the obtainable maps. Therefore a careful evalua- 
tion of the dispersion characteristics of the mixing crystal have to be gauged 
as ultrashort pulses are either investigated or used as a reference. More details 
are discussed in section 4. 




— oo 



S(x,y,Ti,z ) oc (<jAz) 2 \E R \ 2 Io(x,y,Ti,zo) 



(7) 
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2.2 3D field reconstruction 



The method described above is not limited to the intensity reconstruction, 
but can be also implemented for the field reconstruction as can be stated 
from equation 5. This point can also be explained in terms of a holographic 
description of the process, thus bringing to a deeper understanding of the 
imaging reconstruction process. 

As stated in (2^ 13 24), for a plane wave reference the generated wavefront 
through the NL interaction behaves as a conventional hologram recorded at 
frequency u and illuminated with a radiation at 2uo. As a consequence, the 
position, scale, resolution and all the other properties of transformation of the 
reconstructed image can be predicted by means of the ordinary laws of holog- 
raphy. This hologram is recorded and reconstructed at a time, and exists only 
when light propagates inside the crystal. Yet we also point out that we obtain 
the 3D map of our pulse by collecting a set of independent 2D holograms 
by slicing the object pulse at different delays. Nevertheless, according to the 
holographic properties of the SF process, the slicing can be done with the mix- 
ing crystal at any z from the real object to be recovered and the information 
recorded is enough to reconstruct the slice. 

In order to clarify this point, let's assume to perform an experiment in which 
we reconstruct a slice of an object WP (Fig l.a). At any distance z (Fig l.b) 
the reconstructed WP will correspond to the propagated version of the one 
at zq. Let us consider the case when the object WP has a bandwidth small 
enough with respect to the SF bandwidth that our simplified model applies, 
the pulse does not have appreciable angular dispersion 2 and the spatial and 
the temporal evolution can be separated. Under these assumptions, any slice 
recorded in the confocal configuration (see Fig. l.a) could also be reconstructed 
when the mixing crystal is displaced at position z and the field to be converted 
is the propagated one. This is possible by exploiting the holographic properties, 
provided that the detecting system is set to reconstruct the virtual image of 
the slice. This image is placed at a distance 2z far from the mixing crystal 
(see Fig. l.c), its transverse size being identical to that of the object (see (J23h. 
and also note that we are working in the degenerate case, uj\ = c^)- 

The possibility to recover the intensity profile of a virtual image comes from 
the complete wavefront (field) reconstruction arising from the SF process co- 
herence and contained in Eq. 4) 3 . 



2 In the case of angular dispersion, the diffraction drives an effective group velocity 
dispersion also in the vacuum, then the s pat ial evolution cannot be separated from 
the temporal one (see for example (|2lh. ((i^) 

3 We suggest for example the possibility to use an interferometric or heterodyne 
device in order to obtain a complete charaterization of the field. By measuring the 
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Fig. 1. Three reconstruction schemes of a slice obtained from an object WP. In a) 
the slice is reconstructed as a real image (confocal configuration); in b) the slice is 
reconstructed after the object propagated a distance z; in c) the slice corresponding 
to the one imaged in b) is used to reconstruct the virtual image of the one imaged 
in a). 

This important property also allows to get intensity maps with a wider dy- 
namical range. For example, in the cases when the object is so intense so that 
the undepleted pump approximation does not hold, we can get the WP profile 
by displacing the NL crystal to a plane where the intensity is reduced and 
then recover the intensity profile at the desired plane zq by suitably moving 
the imaging system (Fig l.c)). 



3 Experimental results 



We prepared several experiments in order to prove the possibility to recover 3D 
maps of short WPs and to show their holographic properties. The experimental 
set-up is sketched in Fig. 2. 

The 1 ps pulses of a Nd:glass laser source (TWINKLE, Light Conversion, 
wavelength 1055 nm) are splitted on two lines by means of a beam splitter. 
In the first line, the laser pulse is frequency doubled in a KDP crystal and 

field instead of the intensity, a remarkable increment of the dynamical range is also 
possible. 
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UV imaging system Telescope 4x f=+500 mm 

BB ° \ r LBO 




TWINKLE 
(lps, 1055 nm) 



Fig. 2. A sketch of the experimental set-up used to retrieve the space-time intensity 
profiles. 

focused on a 15mm-long litium triborate (LBO) crystal. For pulses about 1 
/xJ in energy , spatial solitons are formed in this crystal because of the optical 
parametric generation process (26). The gaussian temporal profile causes the 
spatial soliton to be formed only in the central part of the pulses where the 
intensity is higher (0)- This leads to a non-trivial space-time structure in 
the output pump wave packet (jH), which we will consider as the object. 
Furthermore, the object has been magnified 4 times by means of a two-lens 
telescope imaging the LBO exit face into the NL mixing crystal ( a 100 fim 
thick, (3 barium borate crystal, BBO). The beam expansion has been necessary 
to: i) reduce the beam intensity and therefore fulfill the undepleted mixing 
requirement; ii) to avoid information loss related to the finite spatial resolution 
of the mixing process. In the second line, a 200-fs reference pulse is produced at 
the wavelength of 527.5 nm by means of a second-harmonic pulse compressor 
([29I ) and expanded 3 times by means of a telescope. Both lines are recombined 
in the thin BBO crystal, cut and oriented to generate the non-collinear SF 
from the object and reference beams. An external incidence angle of ~ 6.5° 
between the propagation directions has been chosen. The delay between the 
two pulses can be varied by means of a suitable delay line placed on the 
reference pulse line. The non-collinear SF radiation is spatially selected by an 
aperture, then selected in frequency by means of coloured filters and the plane 
of the BBO is imaged onto the CCD sensor (PULNIX TM6CN). Finally, both 
the BBO crystal and the imaging system could move independently on a rail 
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Fig. 3. a) Isointensity map of the object pulse as reconstructed from experimental 
data. The rendering has been performed with the Linux-based program GOPEN- 
MOL. b) A contour map of the (x, r) section of the same pulse. 

along the optical axis. 

At first we focused the object inside the BBO crystal and adjusted the position 
of the imaging system as in Fig. l.a). By recording the spatial profile of the 
SF radiation as a function of the reference pulse time lag (steps of At^ = 66.6 
fs has been used, where Ar^ = Ti — Tj_i), we have retrieved the 3D isointensity 
maps of the object pulse at the exit face of the LBO crystal. 

Fig. 3. a shows three different intensity levels of the pulse in the (x, y, r) space, 
while the corresponding contour plot of the (x, r) plane section is depicted 
in Fig. 3.b. The space-time maps clearly show that the pulse is formed by a 
spatially focused structure followed by a diffracted tail. The accuracy of this 
reconstruction has been successfully tested by comparing the experimental 
plots with the results of a 3D-numerical simulation of the formation process 
of the object pulse (see fjiih ). 

In order to check for the validity of the holographic interpretation in this 
process, we have reconstructed the same WP of Figure 3 by moving the NL 
crystal a distance z far from the previous position and by scanning the position 
of the imaging system (CCD+lens) along the SF propagation direction in 
order to find out the position of the (virtual) reconstructed image. First we 
fixed the time lag between the object and the reference and selected a slice 
of the object corresponding to a narrow focused spot of ~ 40/im in diameter 
inside the BBO crystal (corresponding to case a) in Fig. 1). Furthermore, the 
displacements z of the BBO crystal have been chosen large enough to ensure 
that the propagated wavefronts had lost the transverse structure. As a rough 
estimate, for the focused part of the object, the Rayleigh range is about 7 mm 
long, while we spanned distances from -15 mm to 15 mm. In Fig. 4, data show 
the position of the virtual image against the position of the real object from 
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Fig. 4. The positions of the virtual images reconstructed as described in Fig. 1 c) for 
several positions of the real object with respect to the BBO crystal. The two inserts 
show the real image (z = 0, above left) and the virtual image as reconstructed after 
the maximum propagation distance (15 mm, bottom right). All the distances are 
measured with respect to the position of the BBO crystal (positive values, image 
plane positions beyond the crystal in the direction of the beam propagation). The 
points fit a straight line with a slope of about 2, as expected from holography. 

the crystal (see Fig. l.c)), showing a remarkable fit to a straight line which 
slope is close to the value of 2, according to what discussed in the precedent 
section. The error bars shown in the figure indicate the estimated uncertainties 
in the reconstructed image plane position, measured by scanning the whole 
Rayleigh range and by extracting the position of the central point. In the same 
figure the two insets show two isointensity maps recovered for the case when 
the imaging system is focusing directly the BBO plane (z=0) and when the 
BBO crystal has been moved a distance z=15 mm (the maximum propagation 
distance we imposed to the BBO). The agreement between the two intensity 
profiles proves the reliability of our method to work with virtual images in 
recovering 2D slices of WPs like those we used here. 

In order to verify that we have always been operating in the undepleted regime, 
we have measured the dependence of the generated SF field as a function 
of the object-pulse energy at z = 0. The results are presented in Fig. 5, 
where the peak fluence (as registered from the CCD images) is plotted vs 
the transmission of the neutral-density filters that attenuate the object. The 
costant slope confirms the absence of any saturation in the SF process. Data 
also indicate that a slight overestimate of the background has been done during 
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Fig. 5. : The intensity measured by the CCD (arbitrary units) plotted against the 
transmission of the filters used to reduce the object intensity. The linear behaviour 
of Eq. 6 confirms that the measurements were performed in the undepleted-field 
regime. 

the profile acquisition (zero SF field is found for 15 % filter transmission). 



4 Features of the technique 

In the previous sections we have proved that the described technique is able to 
perform the 3D mapping of objects similar to those we used for our measure- 
ments. Anyway, one can expect that in general the reliability of the technique 
to recover the slices could be affected by particular effects dealing with the 
interaction geometry, the mixing crystal, and the features of the interacting 
pulses, in particular for broad-band, ultra-short and chirped pulses. Therefore 
we discuss here three main features of our technique, namely, the duplication 
bandwidth of the mixing crystal, the interaction angle and the reference pulse 
shape, that bring to devise some possible limitations to the fidelity of the 
technique. 

4-1 Space-time resolution 

As we briefly mentioned above, the interaction model we describe does not 
take into account the actual finite bandwidth of the mixing crystal, limiting 
the spatiotemporal resolution of the maps. In fact, in our experiment the BBO 
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crystal has been chosen thin enough that the converted bandwidth was large 
compared to the object one. 

As a matter of fact, the spatio-temporal resolution of the slices is dictated 
by the maximum range of angles and frequencies over which the conversion 
is effective (for very broad-band objects the angular blurring that could arise 
from the phase matching condition should be considered). The key parameter 
here is the maximum phase mismatch between the fundamental and the SF 
waves at which the conversion efficiency vanishes, Ak max . On the basis of the 
existing literature (j3ol). we have estimated the maximum bandwidth converted 
by the crystal as the FWHM of the efficiency curve, yielding to a temporal 
bandwidth of about 630 cm -1 and an angular bandwidth of about 510 cm -1 . 
By considering only the angles for which the conversion efficiency is higher 
than 1/e of the maximum value (jU), we obtain the remarkable resolution of 
details approximately 10 /zm in size (100 lines/mm). We can notice that the 
holographic interpretation of the wavefront reconstruction leads to an easy 
description of the imaging process. Furthermore, our optical system could be 
used in principle to achieve a microscopy of ultrashort objects shorter than 
the one we used here fllih . although the limit of very shot WP is detrimental 
for the holographic reconstruction of the phases. 



4-2 Interaction angle 



A point to be discussed concerns the influence of the noncollinear geometry and 
pulse chirp on the reconstruction of the hologram. In principle the holographic 
interpretation is strictly valid for collinear geometry only (when fc 2 — 2/cx) 
(j23|). Although this condition is not strictly fulfilled in our experiment, we 
have maintained the interaction angle small enough to make this disturbance 
negligible. However, a temporal chirp in the object pulse could give rise to 
a spatial phase distortion of the SF wavefront (|32j). Although this does not 
affect the results obtained with the confocal configuration (see Fig. l.a), it 
could distort the maps obtained in the holographic one (Fig. l.c). We expect 
that this effect is really rampant only when complicated ultrashort pulses with 
strong chirp are considered, or large interaction angle are employed. We are 
confident that our system was operated far from this limit, since, as we checked 
from the data, our experimental holograhic maps do not show any appreciable 
enhancement of the astigmatism. 

As it is well known, the non-collinear SF scheme is largely used to get single- 
shot autocorrelation traces of ultrashort pulses (j33h. In fact, because of the 
geometry of the interaction, it is easy to show that the time of the interaction 
between the object and referencepulses depends on the transverse coordinate 
of the intersecting planes (see (jj, pag. 426-428)). Therefore, if the interact- 
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Fig. 6. Space time map obtained with an external angle of 16 deg between the 
reference and object pulse directions of propagation. The dashed line represents the 
map skewing as calculated from equation 8. 

ing angle is large, we expect that distortions of the spatial profile along this 
direction could affect the space-time maps. However, it turns out that in the 
confocal configuration (i.e. as in Fig. l.a) the effect of the interaction angle is 
merely that of skewing the intensity map, so that the actual time-axis of the 
object pulse is not parallel to that measured on the delay line. This is evident 
in the map depicted in Fig. 6, obtained with an external crossing angle of 



In the paraxial approximation, the slope of the object-pulse time-axis on the 
space-time map is given by: 



where 7 is the angle between the propagation directions of the object and the 
reference wavepackets inside the crystal, and n is the refractive index of the 
medium. It is easy to recognize this formula as the calibration relation of the 
non-collinear, single-shot intensity autocorrelator (jssl : El). 



~ 16°. 
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4-3 Reference pulse shape 



Finally the last point to be discussed concerns the shape of the reference pulse. 
In our experiments it was only slightly affected b y th e presence of satellites 
because of the generation via the pulse compressor 1)291 ) . the amplitude of these 
structures being small enough to be considered negligible to our aim. Actu- 
ally, the recovered intensity profile of the object would be distorted in case 
the reference pulse has a structure more complicated than the single peak en- 
velope, as Eq. 6 clearly points out. This suggests that a careful measurement 
of the reference shape (for example carried out by means of an autocorrela- 
tion technique) can be used to retrieve the real intensity profile by means of 
a deconvolution procedure. We point out, however, that any deconvolution 
unavoidably introduces an extra spurious noise which could degrade the final 
quality of the mapping (0)). 



5 Conclusions 



We have shown that the 3D intensity maps of optical WPs with a com- 
plex space-time structure can be retrived by an optical gating technique. The 
method allows the reconstruction of the WP in its comoving reference frame 
and, by exploiting the holographic properties of a slightly non-collinear de- 
generate, sum-frequency process, we have shown that a complete amplitude 
and phase reconstruction is actually obtained. The maps are obtained by suit- 
ably imaging the second harmonic radiation obtained by cross- correlating an 
object pulse with a much shorter plane wave-packet delayed in time. The 
holographic properties of the generated radiation have been carefully tested 
experimentally, and the distorsions of the intensity maps introduced by the 
non-collinear geometry have been discussed in detail. Finally, theoretical con- 
siderations point out that the ultimate resolution of our set-up is in the order 
of 100 lines/mm in space, and of a few fs in time. However, the choice of the 
reference pulse limits the actual time resolution to about 200 fs. 

We forsee that the developed technique will benefit all the fields where a 
space-time mapping of light pulses is relevant, such as the investigations on 
the resha ping of ultrashort pulses propagating in non-linear materials (flbl 
Hll: lit Eit Il9f). Moreover, the holographic features of the technique might be 
exploited to fully reconstruct the field of an object pulse and its evolution 
during the propagation. 
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